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A two-dimensional system of non-locally coupled complex 
Ginzburg-Landau oscillators is investigated numerically for 
the first time. As already known for the one-dimensional case, 
the system exhibits anomalous spatio-temporal chaos char- 
acterized by power-law spatial correlations. In this chaotic 
regime, the amplitude difference between neighboring ele- 
ments shows temporal noisy on-off intermittency. The system 
is also spatially intermittent in this regime, which is revealed 
by multi-scaling analysis; the amplitude field is multi-afline 
and the difference field is multi-fractal. Correspondingly, the 
probability distribution function of the measure defined for 
each field is strongly non-Gaussian, showing scale-dependent 
deviations in the tails due to intermittency. 



Assemblies of mutually interacting dynamical 
units are ubiquitous in nature. As models for 
them, systems of coupled simple dynamical ele- 
ments, e.g. limit cycle oscillators or chaotic maps, 
have been studied extensively. Recently, a new 
class of coupled systems, namely, a system with 
non-locally coupled dynamical elements, was in- 
troduced and found to exhibit anomalous spatio- 
temporal chaos characterized by power-law spa- 
tial correlations. The preceding studies on this 
system have been done in one dimension, but 
the mechanism for the appearance of this spatio- 
temporal chaos seems to be universal, and it is 
also expected in higher dimensions. In this paper, 
a two-dimensional system of non-locally coupled 
oscillators is investigated for the first time. As 
in the one-dimensional case, the system is found 
to exhibit anomalous spatio-temporal chaos, ac- 
companied by several distinctive features specific 
to this chaotic regime, i.e., power-law spatial cor- 
relations, noisy on-off intermittency, and multi- 
scaling properties. 



I. INTRODUCTION 

Assemblies of dynamical elements coupled with each 
other are widely seen in nature. Simplified models of such 
systems, e.g. coupled limit cycle oscillators or chaotic 
maps, have played important roles not only in modeling 
such systems realistically, but also in understanding the 
varieties of possible behavior of systems far from equilib- 



rium. Many important concepts, such as pattern forma- 
tion or spatio-temporal chaos, have been extracted from 
the detailed studies on such models. 

The interaction between the elements are usually as- 
sumed to be attractive, and of mean-field type in a wide 
sense; each element feels the mean amplitude of its neigh- 
boring elements, and driven by the difference between its 
amplitude and the mean amplitude, in such a way that 
the amplitude differences between the elements decrease. 

The diffusive coupling is a representative limiting case. 
Each element interacts strongly with its nearest neigh- 
bors, so that the amplitude field of the system is always 
continuous and smooth. It is well known that some dif- 
fusively coupled systems of dynamical elements, such as 
the complex Ginzburg-Landau equation, exhibit spatio- 
temporal chaos 

The opposite limiting case is the global coupling, or 
the mean field coupling in the narrow sense. Each ele- 
ment feels the mean field of the entire system, and is thus 
coupled to all the elements with equal strength. The am- 
plitude field becomes statistically spatially homogeneous, 
and the notion of space is lost. It is known that systems 
with global coupling generally show some typical behav- 
ior, e.g. clustering and collective chaos ||,f§. 

In Ref. P], Kuramoto introduced an intermediate 
system between the above two limiting cases, namely 
a system of non-locally coupled elements. Our sub- 
sequent numerical simulations of one-dimensional non- 
locally coupled systems with various elements revealed 
that such systems generally exhibit anomalous spatio- 
temporal chaotic behavior, which cannot be seen in the 
above two limiting cases. In this chaotic regime, the am- 
plitude field becomes fractal, and the spatial correlation 
of the amplitude field shows power-law behavior on small 
scales. Furthermore, the fractal dimension and the ex- 
ponent of the spatial correlation vary continuously with 
the coupling strength. Later, we developed a theory that 
can explain the fractality of the amplitude field and the 
power-law behavior of the spatial correlation based on 
a simple multiplicative stochastic model HQ. Such a 
model is frequently employed in describing the noisy on- 
off intermittency phenomena |S-10 found in many phys- 
ical systems, and this implies that our system should also 
exhibit this type of temporal intermittency. Induced by 
the temporal intermittency, our system is also spatially 
intermittent. In order to study this, we expanded our 
analysis of the amplitude field into more general q-th. 
structure functions, and found that the amplitude field 
shows multi-affinity. We also introduced multi-fractal 
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analysis of the difference field of the original amplitude 
field, inspired by its seeming similarity to the intermit- 
tent energy dissipation field in fluid turbulence [6|]Tf|] . 

All our previous studies have been done in one- 
dimensional systems. However, our previous theory 
does not require the systems to be one-dimensional, and 
spatio-temporal chaos with power-law structure func- 
tions is also expected in higher dimensions. In this pa- 
per, we study a system of non-locally coupled complex 
Ginzburg-Landau oscillators in two dimensions for the 
first time, and investigate its anomalous spatio-temporal 
chaotic regime. Some attention is paid to the multi- 
scaling properties of the intermittent amplitude and dif- 
ference fields. 
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where d is the dimension of the space. When the sys- 
tem is isotropic, the kernel G becomes a function of the 
distance r := \r' — r\, and is expressed as 



G(r) oc exp(— j\r\) d=l, 
K (rf\r\) d=2, 
exp(— 7|r|) 
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where Kq is the modified Bessel function. The constant 7 
gives the inverse of the coupling length, and is calculated 

as 



II. MODEL 

As proposed by Kuramoto ||, the non-local coupling 
naturally appears in the following plausible situation. 
Consider an assembly of spatially distributed dynamical 
elements, e.g. cells. Each element is assumed to interact 
indirectly with other elements through some (e.g. chem- 
ical) substance, which diffuses and decays much faster 
than the dynamics of the element. Such a situation will 
be described by the following set of equations: 

X(r,t)=F(X(r,t))+K-A(r,t), (1) 
eA(r, t) = -T)A(r, t) + DV 2 A(r, t) + X(r, t), (2) 

where X is the amplitude of the element, F is the dy- 
namics of the amplitude, and A is the concentration of 
the substance with decay rate 77 and diffusion rate D. 
The substance A is generated at a rate proportional to 
the amplitude X, and the amplitude X is affected by the 
substance A with a coupling matrix K. The parameter 
e determines the ratio of time scale of the elements to 
that of the substance, and is assumed to be very small. 
Namely, the dynamics of the substance is much faster 
than that of the elements. 

Now, let us consider the e — + limit and eliminate the 
dynamics of A adiabatically. Putting the left-hand side 
of Eq.(||) to 0, we can solve the equation for A as 



(10) 



A(r,t)= f dr'G(r' -r)X(r',t), 



(3) 



where G(r' — r) is a kernel that satisfies 

(77 - DV 2 )G{r' - r) = S(r'). (4) 

By inserting Eq.(||) into Eq. ([!]), we obtain the following 
system of non-locally coupled dynamical elements: 



X(r, t) = F(X(r, t)) + K • / dr'G{r' - r)X(r', t). (5) 
The kernel G(r' — r) can be solved as 



Each G(r) must satisfy the normalization condition 
J G(r)d d r = 1. Since we treat a two-dimensional sys- 
tem, we use Eq.(||) for G(r) hereafter. 

As elements, we use complex Ginzburg-Landau oscil- 
lators. They are the simplest limit cycle oscillators that 
can be derived by the center-manifold reduction tech- 
nique from generic oscillators near their Hopf bifurcation 
points jjj. The corresponding non-locally coupled sys- 
tem is given by the following equation for the complex 
amplitude W: 

W(r, t) = W-(l + ic 2 )\W\ 2 W + K(l + ia){W - W), 

(11) 

where K is the coupling strength, C\,C2 are real param- 
eters, and the non-local mean field W is given by 



W{r,t)= dr'G{r' -r)W(r',t). 



(12) 



This is the non-local complex Ginzburg-Landau equation 
introduced by Kuramoto || as the first concrete example 
of non-locally coupled systems. 



III. ANOMALOUS SPATIO-TEMPORAL CHAOS 

In the numerical simulations presented here, we as- 
sume the total system to be a square lattice with both 
sides unit length long. The elements are placed on the 
lattice sites, and periodic boundary conditions are as- 
sumed. We use A^ 2 = 512 2 ~ 1024 2 elements, and fix 
the coupling length 7 _1 at 1/8. The non-local mean field 
is easily calculated by using the FFT technique, since it 
is simply a convolution of the amplitude field with the 
kernel (||). We fix the parameters c\ and C2 at —2 and 2 
respectively. These are the standard values already used 
in our previous one-dimensional simulations. 



2 



In Fig.[l]-[| typical snapshots of the real part X(x,y) 
of the complex variable W(x,y) are shown for three dif- 
ferent values of coupling strength K. Since we obtain 
similar figures for the imaginary part Y(x, y) by symme- 
try, we use X(x,y) in the following analysis and call it 
amplitude field. The amplitude field at K = 1.05 is con- 
tinuous and smooth, while at K = 0.65 it seems to be dis- 
continuous and disordered, although not completely ran- 
dom. The amplitude field at the intermediate coupling 
strength K — 0.85 looks somewhat more complex and in- 
triguing; it is composed of intricately convoluted smooth 
and disordered patches of various length scales. This is 
the anomalous spatio-temporal chaotic regime that our 
interest is focused on. 




FIG. 1. Snapshot of the amplitude field X(x,y) at 
K — 1.05. Darker dots indicates larger amplitudes. 



IV. SPATIAL CORRELATION FUNCTION 

Let us examine the spatial correlation function 
first. Figures ||(a)-(c) show spatial correlation functions 
C(x,y) := (X(0, 0)X(x, w)) corresponding to the ampli- 
tude fields shown in FigsJ^-^. Each correlation function is 
clearly rotationally symmetric due to the isotropy of the 
system. As the amplitude field becomes disordered, the 
correlation function becomes steep, and the center of the 
graph, which corresponds to the self-correlation C(0, 0), 
becomes peaked. 

In Ref. ||, the anomalous spatio-temporal chaotic 
regime was characterized by the power-law behavior of 
the spatial correlation function in small distance: 

C(l) := (X(0)X(/)) ~C - C x r (I « 1), (13) 




FIG. 2. Snapshot of the amplitude field X(x,y) at 
K = 0.85. 




FIG. 3. Snapshot of the amplitude field X(x,y) at 
K = 0.65. 
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where C$,C\ are constants and a is a non-integer 
parameter-dependent exponent. 

To confirm if this power-law behavior also holds in two 
dimension, we calculated the radial correlation function 
C(l) = (X(r)X(r + l)) (\l\ = I) along a straight line in a 
certain direction)^, and estimated the best fitting param- 
eters Co and C\. Figure | shows ln^ vs. In [C - C{1)] 
for several values of the coupling strength K. For each 
coupling strength, the experimental data are almost on a 
straight line, and the power-law behavior is evident. The 
exponent a of the power law varies continuously with 
the coupling strength. Although not shown in the figure, 
the correlation function C(l) is continuous at the origin 
I = for K > 0.85, but discontinuous at K = 0.80. 
There appears a finite gap between the self-correlation 
C(0) and the correlation between the nearest-neighbor 
elements limi^+o C(/) ~ Cq. This means that the in- 
dividual motion of the element becomes so violent that 
the amplitude field is no longer continuous statistically. 
In Ref. fjj], this transition point was identified with the 
blowout bifurcation point in the on-off intermittent dy- 
namics of the amplitude difference between nearby ele- 
ments. 

Thus, the anomalous spatio-temporal chaos in two 
dimension is also characterized by power-law behavior 
of the spatial correlation function with a parameter- 
dependent exponent. 



V. NOISY ON-OFF INTERMITTENCY 

In Ref. |,|, we clarified that the scaling behavior 
of the spatial correlation is a consequence of underly- 
ing multiplicative processes of amplitude differences be- 
tween neighboring elements. We described this process 
by a multiplicative stochastic model, and related the ex- 
ponent a of the spatial correlation with the fluctuation of 
the finite-time Lyapunov exponent of the element. Such 
a model is essentially identical with those used in de- 
scribing noisy on-off intermittency |^— [Tc|] , which indicates 
that our system also exhibits this type of temporal inter- 
mittency. Actually, the finite-time Lyapunov exponent 
of the complex Ginzburg-Landau oscillator can fluctuate 
between positive and negative values, and neighboring os- 
cillators are subjected to only slightly different non-local 
mean field. Therefore, the conditions for the appearance 
of noisy on-off intermittency are satisfied in our system. 

Now, let us confirm this in our system numerically. 
The coupling strength is set at K — 0.85. Figure |6| shows 
typical time sequences of amplitude differences AXi(t) 
and AX2(t). The distance between the elements is 512 -1 
for AXi(t), and 64 _1 for AX 2 (t), respectively. Strong 



1 We mainly used the (0, 1) or the (1, 
are independent of the direction. 
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FIG. 4. Short-range spatial correlation functions C(x, y) at 
(a) K = 1.05, (b) K = 0.85, and (c) K = 0.65. 
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FIG. 5. Power-law behavior of the correlation functions; 
log-log plot of Go — C(Z) vs. I obtained for several values of 
the coupling strength K. 
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intermittency of the signals is apparent. It can be seen 
that AX2(t) shows more frequent bursts than AXi(t), 
reflecting that AX 2 (t) is subjected to larger fluctuations 
than AXi(t). 

We can confirm that these intermittent signals are ac- 
tually noisy on-off intermittent by calculating the lam- 
inar length distribution. The laminar phase is defined 
as a successive duration where the absolute value of the 
difference does not exceeds a certain threshold. Here 
we choose 0.5 as the threshold value. Figure shows 
laminar length distributions R(t) obtained from AXx(t) 
and AX2(t). The characteristic shape of the distribu- 
tion R(t), i.e., the power-law dependence on t with slope 
—3/2 for small t, and the exponential shoulder seen in 
the large t region, clearly indicates that the signals are 
actually noisy on-off intermittent. The shoulder reflects 
broken scale invariance due to the additive noise. As ex- 
pected, the shoulder of appears at a smaller value 
of t than that of AX^t) §|o|. 
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1000 2000 3000 4000 5000 6000 
t 

FIG. 6. Typical temporal evolution of the amplitude dif- 
ferences. The distance between the elements is 512" 1 for 
AXi(t), and 64" 1 for AX 2 {i). 



VI. MULTI-SCALING ANALYSIS 

The notion of multi-scaling, i.e., multi-affinity and 
multi-fractality, have been employed successfully in char- 
acterizing complex spatio-temporal behavior of various 
phenomena, such as velocity and energy dissipation fields 
in fluid turbulence rough interfaces in fractal sur- 

face growth fl3]] , nematic fluid electro-convective turbu- 
lence 14 1 , financial data of currency exchange rates (jl5j , 
and even in natural images Jl6| . In Ref . ||0 , we intro- 
duced multi-scaling analysis to our system for the one- 
dimensional case, inspired by the seeming similarity of 
the amplitude and difference fields in our system to the 
velocity and energy dissipation fields in fluid turbulence. 
Here, we attempt the multi-scaling analysis for the two- 
dimensional case. 




FIG. 7. Distributions of the laminar length obtained from 
the time sequences AXi (t) and AX2 (t) shown in Fig.^| 



First, we introduce the difference field Z{r) as 



Z(r) := \WX(r)\ 



:f) 2 +(f)' <»> 



which emphasizes the edges of the original amplitude field 
X(r) 2 . This quantity is an analogue of the energy dissi- 
pation field in fluid turbulence. Figure |8| shows a typical 
snapshot of the difference field Z(x, y) at K = 0.85, cor- 
responding to the amplitude field shown in FigJ|. The 
intermittency underlying the original amplitude field is 
now apparent. 

We then introduce the following quantities as measures 
for the amplitude field X(r) and the difference field Z(r): 



h(r;l) := \X (r + I) - X (r)\ 

m(r;l) := [ Z(r')d 2 r', 
Jsfra) 



(15) 
(16) 



where \l\ = I, and the domain of integration S(r;l) is 
a square of size I placed at r. The first quantity is a 
difference of the amplitude field X(r) between two points 
separated by a distance of I, and the second quantity is 
a volume enclosed by the difference field Z(r) and the 
square S(r;l). 

According to the multi-fractal formalism, two types of 
partition functions are defined as 



z q h (i) ■■- 



(h(l) g ) 



1 



M(l) 



(17) 



2 Here, the differential should not be interpreted literally. 
We always use finite difference in the actual calculation, e.g. 
(X(x + Ax, y) — X(x, y))/ Ax with sufficiently small Ax, and 
this is important to observe the multi-scaling behavior [hll. 



5 




FIG. 8. Snapshot of the difference field Z(x,y) at 
K — 0.85, corresponding to the amplitude field shown in 
Fig|. 



N(l) 

^(0:^(0M0 9 HE m ( r ^) 9 , (is) 

j=l 

where Z?(l) is calculated along a certain straight line 
in some direction as in the case of the previous spatial 
correlation function, while Z q n (V) is calculated over the 
whole system. r j is either the position of the line segment 
or the position of the square. M(l) is the number of 
line segments of length I that are needed to cover the 
entire line, and N(l) is the number of squares of size I 
that are needed to cover the whole system. The function 
Zt(l) is called structure function in the context of fluid 
turbulence §,|||. 

When the measures have scaling properties, the parti- 
tion functions are expected to scale with I as Zi(t) ~ l^ q ~) 
and Z^il) ~ l T ^ q \ Furthermore, if these exponents ((q) 
and r(q) depend nonlinearly on q, the corresponding 
measures h(l) and m(l) are called multi-affine and multi- 
fractal, respectively H jl^ijlj . 

For the one-dimensional case, we already know that 
the amplitude field is multi-affine and the difference field 
is multi-fractal. Moreover, our previous theory predicts 
the following form for the scaling exponent of the 
amplitude field: 

C(?) = <7 (0 < q < /?), P(P<q), (19) 

where (5 is a positive constant determined by the fluctua- 
tion of the finite-time Lyapunov exponent of the element, 
and is related to the slope of the probability distribution 
of h(l) pljLl]]. This is the simplest form of multi-affinity, 



and sometimes called bi-fractality |l2l]n. The same form 
of the scaling exponent ((q) is also expected in two di- 
mension, since our previous theory imposed no restriction 
on the dimensionality of the system. 

For the scaling exponent r(q) of the difference field, 
we have not been able to develop a satisfactory theory 
yet. Numerical results in one-dimensional systems sug- 
gest that r(q) also depends nonlinearly on q, and the 
difference field is multi-fractal with a rather simple func- 
tional form for r{q) However, the scaling exponent 
r(g) for the two-dimensional system may be different 
from that for the one-dimensional case, since Z^il) is 
defined depending on the dimensionality of the system, 
while Z?(l) is always defined along a one-dimensional 
line. 
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FIG. 9. Partition functions Z^(l) vs. I for several values of 
q. The top line corresponds to q = 0, and the bottom one to 
q = 5, at intervals of 0.5. 



Let us proceed to the numerical results now. The cou- 
pling strength is fixed at K — 0.85 hereafter, where the 
system is fully in the anomalous spatio-temporal chaotic 
regime. 

Figure || shows the partition function Z^(l) obtained 
for several values of q. Each curve depends on I in a 
power-law manner for small I, and its exponent increases 
with q. The dependence of the scaling exponent £(<?) on 
q is shown in Fig|ll]. The £(q) curve is a strongly nonlin- 
ear function of q, and the multi-scaling property of the 
amplitude field is evident. Furthermore, the £(q) curve 
has a bi-linear form as expected in Eq.(fl9|), although a 
sharp transition is absent due to the limited number of 
oscillators. From the large q behavior of the exponent, 
we can roughly estimate the value of (3 as ~ 0.48. 

Thus, the amplitude field turns out to be multi-affine, 



3 There is some confusion in terminology due to historical 
reasons. A word 'bi-affinity' will be more appropriate if it 
exists. 
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FIG. 10. Partition functions Z^(l) vs. 
of q. The bottom line corresponds to q = 
to q — 5, at intervals of 0.5. 
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I for several values 
0, and the top one 



and the behavior of the scaling exponent is the same as 
that for the one-dimensional case. 

Figure |l0| shows the partition function Z^(l) for sev- 
eral values of q. It is clear that each curve shows a power- 
law dependence on I. The width of the region where the 
power law holds seems much wider than the previous 
case. The scaling exponent r(q) is plotted with regard 
to q in Fig.[l2| The corresponding generalized dimension 
D(q) :— r(q)/(q — 1) is also shown in the inset. The 
r(q) curve is again a nonlinear function of q, but its de- 
pendence on q does not seem to be so simple as that of 
the <^(q) curve. However, as we conjectured numerically 
for the one-dimensional case [Tl| l , asymptotic linearity of 
the r(q) curve seems to hold. Correspondingly, the D(q) 
curve seems to saturate to a horizontal line D(q) — D(oo) 
quickly. But we can not observe a clear transition to the 
horizontal line as in the previous one-dimensional case, 
which may be due to the limited number of oscillators, 
or the two-dimensionality of the system. 

Thus, the difference field also turns out to be multi- 
fractal. The behavior of the scaling exponent somewhat 
resembles that in the one-dimensional case, but is not 
completely the same. 
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11. Scaling exponent £(q) vs. q. The theoretical curve 
with (3 = 0.48 is compared with the experimental 
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VII. PROBABILITY DISTRIBUTIONS OF THE 
MEASURES 

The multi-scaling properties of amplitude and differ- 
ence fields are consequences of intermittency underlying 
the system. In order to analyze this intermittency in 
more detail, we study here the probability distribution 
functions (PDF) of both measures at each length scale. 

Let us consider the PDFs of the measures h(l) and 
m(Z). It is convenient to use rescaled measures h r (l) :— 
h(l)/l, m r (l) :— m(l)/l 2 and corresponding rescaled 
PDFs P r (h r ; I), Q r (m r \ I). With this rescaling, the peaks 
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q 

FIG. 12. Scaling exponent r(q) vs. q. The inset shows the 
corresponding generalized dimension D(q) vs. q. 
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and the widths of the PDFs become relatively close. We 
show our results in these rescaled variables. 

The scaling exponents ((q) and r(q) are fully deter- 
mined by the dependence of these PDFs on the scale I, In 
Rcf . , we approximated the tails of the PDFs by cer- 
tain functional forms and extracted the scaling exponents 
asymptotically in the small I limit. Here we present the 
numerical results only briefly for the purpose of empha- 
sizing the intermittency of the amplitude and difference 
fields. 

Figure |l3| shows the rescaled PDF P r (h r ;V) of the 
rescaled amplitude difference h r (l) for several values of 
I in log-log scales. Each PDF has a characteristic trun- 
cated Levy-like shape, as we already obtained previously 
for the one-dimensional case; it is composed of a con- 
stant region near the origin, a power-law decay in the 
middle, and a sharp cutoff. Each curve roughly collapses 
to a scale-invariant curve in the constant and power-law 
regions, while the cut-off of the tail moves to the right 
with decreasing of I due to the intermittency. More pre- 
cisely, the cut-off position of the PDF (defined in some 
suitable way) is proportional to I . This gives rise to 
the bi-fractal behavior of the curve, see Ref. for 
detail. 

Our previous theory predicts that the slope of the 
power- law decay is given by —1 — (3 with the constant 
(3 from Eq. (|l9|). This is confirmed in Fig.|l3|, where the 
slope of the power-law decay can be read off as —1.4, 
roughly in agreement with the previously obtained value 
/3 ~ 0.48 from the scaling exponent ((q). The inset 
of Fig. [l3| shows the PDF of the amplitude difference 
h(l) (without taking the absolute value) rescaled by the 
standard deviation in linear-log scales, in order to fur- 
ther emphasize the intermittency of the amplitude field. 
The PDF evolves from nearly Gaussian into intermittent 
power-law with the decrease of scale I. 

Figure [li] shows the PDFs Q r (m r ;l) of the rescaled 
volume m r (l). Their shapes are not so simple as those 
for P r (h r ;l). They are also qualitatively different from 
the PDF for the difference field in the previous one- 
dimensional case , which is responsible for the slightly 
different behavior of the scaling exponent r{q) from the 
one-dimensional case. But we can still see that both tails 
extend, and the distribution widens with the decrease of 
scale I. Namely, as we decrease the observation scale 
I, largely deviated events appear more frequently. It is 
obvious that this intermittency effect gives rise to the 
nonlinearity of the scaling exponent r(g), although the 
precise functional form is difficult to obtain. The inset 
shows the PDF of the measure m(l) rescaled by the stan- 
dard deviation in linear-log scales for several values of 
The PDF gradually gets steeper with the decrease of I 
due to the intermittency. 

Thus, the PDFs of the measures reveal the intermit- 
tency in our system clearly. Especially, the PDF for the 
amplitude difference h(l) has the same shape as that al- 
ready obtained in the one-dimensional case. 




h T 

FIG. 13. Rescaled PDFs P r {h r ;l) for I = \/2, 2a/2, 4^, 
8V2, 16V2, and 32\/2 (x 1024" 1 ). The curve with the slowest 
cutoff corresponds to I — \/2 x 1024" 1 , and the leftmost curve 
with the fastest cutoff corresponds to I = 32\/2 x 1024 -1 . The 
inset shows PDFs of the amplitude difference h(l) (without 
taking the absolute value) rescaled by the standard deviation 
in log-linear scales for I — \fl (the steepest curve) , 4\/2, 16\/2, 
64\/2, and 128\/2 (the nearly quadratic curve) (X1024- 1 ). 




FIG. 14. Rescaled PDFs Q r (m r ;l) for I = 2 (the widest 
distribution with the steepest power-law decay), 8, 32, 128 
(the narrowest) (xl024 _1 ). The inset shows PDFs of m(l) 
rescaled by the standard deviation in log-linear scales for 
the same set of I values. The steepest curve corresponds to 
U 2 x 1024" \ and the curve with nearly quadratic peak 
corresponds to I = 128 x 1024" 1 . 
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VIII. CONCLUSION 

We numerically analyzed a two-dimensional system of 
non-locally coupled complex Ginzburg-Landau oscilla- 
tors. As in the previous one-dimensional case, we found 
an anomalous spatio-temporally chaotic regime charac- 
terized by a power-law behavior of the spatial correla- 
tion function. As expected from our previous theory, 
the amplitude difference between neighboring elements 
exhibits noisy on-off intermittency, giving a microscopic 
dynamical origin for the power-law spatial correlations. 
We performed multi-scaling analysis in that regime, and 
found that the amplitude and difference fields are indeed 
multi-affine and multi-fractal, indicating strong intermit- 
tency underlying the system. By studying the PDFs 
of the measures at each length scale, the intermittency 
was clearly observed as scale-dependent deviations of the 
PDFs in their tails. 

Multi-scaling properties are also known in various phe- 
nomena such as turbulence or fractal surface growth. The 
appearance of similar multi-scaling properties in many 
different systems implies some underlying common sta- 
tistical law leading to such behavior. Further study of the 
intermittency in our system will give more insights and 
hints for the understanding of the multi-scaling proper- 
ties observed in complex dissipative systems. 
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